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ABSTRACT 

We present the galaxy cluster autocorrelation function of 277 galaxy cluster candidates with 0.25 < z < 1.5 
in a 7 deg^ area of the IRAC Shallow Cluster Survey. We find strong clustering throughout our galaxy cluster 
sample, as expected for these massive structures. Specifically, at (z) = 0.5 we find a correlation length of 
ro = 17.40^3 jQ /z~' Mpc, in excellent agreement with the Las Campanas Distant Cluster Survey, the only other 
non-local measurement. At higher redshift, (z) = 1, we find that strong clustering persists, with a correlation 
length of ro = 19. 14^^ ^^ /i"' Mpc. A comparison with high resolution cosmological simulations indicates these 
are clusters with halo masses of ~ IO'^^Mq, a result supported by estimates of dynamical mass for a subset 
of the sample. In a stable clustering picture, these clusters will evolve into massive (lO'^M©) clusters by the 
present day. 
Subject headings: galaxies: clusters: general — cosmology: observations — large-scale structure of the uni- 
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1. INTRODUCTION 

The clustering amplitude of massive galaxy clusters, and 
in particular its dependence on richness or cluster mass, is a 
strong function of the underlying cosmology. Recent theo- 
retical work (e.g., Majumdar and Mohr 2003; Younger et al. 
2006) has demonstrated how cluster surveys can be "self- 
calibrated," providing precise simultaneous constraints on 
both the cosmology and cluster evolution models. Given red- 
shift and approximate mass information for each cluster, re- 
liable cosmological parameter estimation is feasible even in 
the presence of significant, and potentially unknown, evolu- 
tion in cluster physical parameters (e.g., Gladders et al. 2007). 
The ultimate goal of these cosmological studies, a measure- 
ment of the equation of state of dark energy with an accuracy 
competitive with SNe la methods, also requires knowledge of 
the cluster power spectrum or autocorrelation function (Ma- 
jumdar and Mohr 2004; Wang et al. 2004). Such extensive 
sample characterization naturally emerges from mid-infrared 
selected photometric redshift cluster surveys like the IRAC 
Shallow Cluster Survey (ISCS; Eisenhardt et al. 2007, here- 
after E07). 

Previous studies (e.g., Bahcall et al. 2003, and references 
therein) have produced largely local (z < 0.3) galaxy cluster 
clustering measurements with limited baselines for evolution- 
ary studies. The ISCS provides the opportunity to address 
these issues in one of the largest statistical samples of high 
redshift clusters to date. In this Letter we present the first 
measurement of the galaxy cluster autocorrelation function 
extending over more than half the cosmic age of the Uni- 
verse. Such measurements also allow evolutionary connec- 
tions between galaxy clusters and high redshift, highly clus- 
tered galaxy populations to be explored. 

We use a concordance cosmology throughout, with Hm = 
0.3, r^A = 0.7, and Hq = 70 km s"' Mpc"'. For consistency 
with previous studies we report distances, including correla- 
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tion lengths, in units of comoving h ' Mpc, with //q = lOO/i 
kms"' Mpc"'. 

2. IRAC SHALLOW CLUSTER SURVEY 

The ISCS is a sample of 335 galaxy clusters spanning 
0.1 < z < 2 in the IRAC Shallow Survey (ISS, Eisenhardt 
et al. 2004) area of the NOAO Deep-Wide Field Survey (ND- 
WFS, Jannuzi and Dey 1999; Jannuzi et al. in prep) in Bootes. 
The clusters are selected using a wavelet detection algorithm 
which identifies peaks in cluster probability density maps con- 
structed from accurate photometric redshift probability func- 
tions for 175,431 galaxies brighter than 13.3/iJy at 4.5/im 
in a 7.25 deg^ region (Brodwin et al. 2006, hereafter B06; 
E07). The galaxy photometric redshifts, which are key to the 
cluster-finding algorithm, are derived from the joint ISS, ND- 
WFS (DR3) and FLAMEX (Elston et al. 2006) data sets. 

The AGES survey (Kochanek et al. in prep) in Bootes pro- 
vides spectroscopic confirmation for dozens of clusters at 
z < 0.5 (E07). At higher redshift, a multi-year Keck spec- 
troscopic campaign has to date confirmed 10 z > 1 clusters 
(Stanford et al. 2005; Elston et al. 2006; B06; E07). 

Extensive masking, described in B06, is implemented to re- 
ject areas suffering from remnant cosmetic artifacts or data 
quality issues which could affect the robustness of the clus- 
tering measurements. The final unmasked area is 7.00 deg^ 
and contains 320 of the 335 galaxy clusters in the full ISCS, 
of which 277 are within the redshift range (0.25 < z < 1.5) 
considered in this Letter 

3. GALAXY CLUSTER AUTOCORRELATION MEASUREMENTS 

3.1. Angular Correlation Function in Redshift Slices 

The redshift distribution of galaxy clusters is presented in 
Figure 1. A fit of the formA^(z) oc (z/z())"exp(-(z/zo)''), with 
a= 1.5, /3 = 1.7 and zo = 0.7, is overplotted. Cosmological 
studies with this observed redshift distribution must await a 
detailed characterization of the survey selection function, to 
be presented in a forthcoming paper. Galaxy clusters are split 
into two similarly populated redshift bins, 0.25 < z < 0.75 
(136 clusters) and 0.75 <z< 1.5 (141 clusters), using the best 
available redshift information (i.e., spectroscopic redshifts 
where available, photometric redshifts otherwise). These bins 
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span the range where the accuracy, a, = 0.06(1 +z), and reH- 
ability of the galaxy photometric redshifts have been demon- 
strated (B06; see also Brown et al. 2007). 
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Fig. 1. — Observed redshift distribution of galaxy clusters in the ISCS. The 
dashed lines illustrate the two redshift bins considered here, 0.25 < z < 0.75 
and 0.75 < z < 1.5, containing 277 clusters between them. The curve is a fit 
to the distribution and is used in the Limber deprojection. This distribution 
should not be used for quantifying number density evolution without careful 
inclusion of selection biases. 



The angular correlation function (ACF) is parametrized 
here as a simple power law. 



oj(e)=A, 



(1) 



This can be deprojected (Limber 1954) to yield a measure- 
ment of the real-space correlation length, r()(z), over the red- 
shift range spanned by the 2-D sample. 



r2{Ze«)=A^ 



HqH^ iyHz)[x(z)V-^E(z)dz 
c [C;N(z)dz]^ 



where 



_J;^zNHz)[x(z)]'-^E(z)dz 
^=ff- rNHz)[x(z)]'-'<E(z)dz' 



(2) 



(3) 



-f=l + 5,H^ = r(l/2) r[(7- l)/2]/r(7/2), N(z) is the red- 
shift distribution, and E{z) and x(z) describe the evolution of 
the Hubble parameter and the comoving radial distance, re- 
spectively (e.g., Hogg 1999). Since the cluster redshift uncer- 
tainties are much smaller than the width of our redshift bins 
(E07), they have little impact on the redshift distribution and 
the modeling of the spatial correlation function. 

3.2. Results 

We calculate the ACF using both the Landy and Szalay 
(1993) and Hamilton (1993) estimators using 500,000 ran- 
doms to ensure a robust Monte Carlo integration. The re- 
sults are nearly identical with both estimators and we report 
the results obtained with the latter Two independent fitting 
techniques are applied to the data. The standard frequentist 
(or classical) approach is used to simultaneously fit the slope, 
and, through the use of the relativistic Limber equation, the 



correlation length, ro. In calculating these correlation lengths 
we have adopted the N(z) parametrization shown in Figure 1 . 
A Bayesian technique is also employed to directly determine 
the correlation lengths, marginalizing over the slope, 5, sub- 
ject to the weak prior that it is in the range 0.2 < S < 1.8. 
This method is desirable since, despite the large number of 
clusters, it is difficult to simultaneously constrain both ampli- 
tude and slope. Figure 2 shows both the frequentist fits and 
the Bayesian likelihood functions in ro. 
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Fig. 2. — Angular correlation functions in redshift bins at 0.25 < z < 0.75 
(top) and 0.75 < z < 1.5 (bottom). The error bars are estimated via bootstrap 
resampling. The red solid lines show the best-fit ACFs, with con'esponding 
slopes listed in Table 1 . The blue dashed fines show the fits with a fixed slope 
(<5 = 1.0). Marginalizing over the slope produces the likelihood functions in 
ro presented in the insets. The 68% and 90% confidence intervals are shown. 

In this highly clustered population, sample (or cosmic) vari- 
ance dominates over simple Poisson errors. We therefore cal- 
culate the error bars using 100 bootstrap resamplings with re- 
placement. To test for possible systematics across the field 
we divided the field into halves, once north-south and once 
east-west, and for each subfield we computed the ACF. The 
results for all half-fields agree within 1 a, indicating that we 
are not adversely affected by an unidentified bias in the spatial 
selection function. 

Bootstrap simulations of the cluster detection process indi- 
cate that the spurious fraction is less than 10% at all redshifts 
(E07), and spectroscopic observations indicate it is likely 
much lower. Conservatively assuming 5-10% of the sample is 
indeed spurious, and that these are uncorrelated, then at most 
we are underestimating the clustering by « 1 1-23%. 

In this work we adopt the ro values from the Bayesian fits. 



Galaxy Cluster Correlation Function to z ~ 1.5 



though the fits from both methods, presented in Table 1, are 
completely consistent. The correlation amplitude in our low 
redshift bin at z « 0.5 is in excellent agreement with the only 
other measurement at this redshift (Gonzalez et al. 2002). Our 
high redshift measurement, at z « 1, is the first to probe struc- 
ture on the largest scales in the first half of Universe. 

The space densities for these samples and the mean inter- 
cluster distances, dc, are also presented in Table 1. The re- 
lationship between clustering amplitude and d^ predicted in 
a concordance cosmology, and observed in practice (Bahcall 
et al. 2003, and references therein), is only weakly dependent 
on redshift. As shown in Figure 3, the ISCS samples are quite 
consistent with the LCDM predictions between < z < 1.5 
(hashed region) from Younger et al. (2005). 
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Fig. 3. — Plot of ro vs. dc for the present sample, along with several 
measurements taken from the compilation of Bahcall et al. (2003). Our 
results (filled blue circles) are consistent with LCDM predictions between 
< z < 1.5 (hashed region), with the parametrization of Younger et al. (2005, 
dashed line), as well as with the most distant measurement prior to the present 
one, the LCDCS at z ~ 0.5 (open blue circles). Data are from Abadi et al. 
(1998); Bahcall and Soneira (1983); Collins et al. (2000); Croft et al. (1997); 
Gonzalez et al. (2002); Lee and Park (1999); Peacock and West (1992). 



4. DISCUSSION 

A key theoretically predictable cluster observable is the cor- 
relation function as a function of halo mass. In simulations 
the halo mass, Mzoo, is defined as the mass inside the radius 
at which the mean overdensity is 200 times the critical den- 
sity. We compare our clustering results with the Younger et al. 
(2005) analysis of the Hopkins et al. (2005) high-resolution 
cosmological simulation, which had a 1500 /i~' Mpc box 
length, an individual particle mass of 1.8 x IO^Mq, and a 
power spectrum normalization of a^ = 0.84. We infer that 
the ISCS cluster sample has average log[M2oo/A^o] masses of 



13.9!!jjand' 



13.8!{J ;f at Zeff = 0.53 and 0.97, respectively. 



Direct dynamical masses for the 10 z > 1 clusters presented 
in £07 yield a largely consistent distribution of masses (Gon- 
zalez et al. in prep). 

The observed constancy of ro for massive galaxy clusters 
out to z = 1 is a robust confirmation of a key prediction from 
numerical simulations, reflecting the relative constancy of the 
mass hierarchy of clusters with redshift (Younger et al. 2005). 
That is, the N most massive clusters at one epoch roughly 
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Fig. 4. — Comoving coirelation lengths for the ISCS clusters, along with 
other cluster samples at lower redshift (offset slightly in z for clarity) and 
highly clustered galaxy samples at higher redshifts (see the text for a full 
description and references). Two evolutionary models are overplotted, in- 
cluding the Fry (1996) biased structure formation model (shaded region and 
dotted lines for Zgff = 0.53 and 0.97, respectively), and a simple stable clus- 
tering model (dashed lines). 



correspond to the N most massive clusters at a later epoch, 
and therefore have similar clustering. 

In Figure 4 we plot a compilation of recent clustering am- 
plitudes for various cluster surveys, as well as for highly 
clustered galaxy populations, including FRII radio galaxies 
(Overzier et al. 2003), EROs (Brown et al. 2005; Daddi 
et al. 2001, 2004), ULIRGs (FaiTah et al. 2006; MagHoc- 
chetti et al. 2007), and SMGs (Blain et al. 2004). Since 
optically-selected QSOs are more modestly clustered (Coil 
et al. 2007; Croom et al. 2005; Myers et al. 2006; Porciani 
et al. 2004), and therefore reside in considerably less massive 
halos {^ 10'^- IO'^Mq) than the ISCS clusters, they are not 
included in Figure 4. 

Following Moustakas and Somerville (2002), we overplot 
the halo conserving model of Fry (1996) normalized to our 
two measurements in order to explore possible evolutionary 
connections with structures at other redshifts. The shaded 
area (dotted lines) shows the 1 a region for the Zeff = 0.53 
(0.97) measurement. In this model, representative of a class 
of merger-free biased structure formation models, the ISCS 
clusters will evolve into typical present-day massive clusters, 
such as those in the SDSS, APM or Abell surveys. In the 
stable-clustering picture, in which clustering is fixed in phys- 
ical coordinates (Groth and Peebles 1977, dashed lines), the 
Zeff = 0.97 ISCS clusters grow into the most massive clusters 
in the local Universe, typically identified in X-ray surveys. 

Most of the plotted high redshift galaxy clustering measure- 
ments are rather uncertain due to both small number statis- 
tics and poorly known redshift distributions. Mindful of this 
caveat, we observe that FRII radio galaxies, some ERO sam- 
ples, and z '^ 2 ULIRGs have clustering consistent with that 
seen in the ISCS clusters in either model. Clearly these pop- 
ulations trace very massive halos (> IQ^^Mq). As a cau- 
tion against overinterpretation, however, we note that only 
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Table 1 . Amplitude of Correlations for ISCS Clusters 



Az Zg,ff N logA^" 



'■o 
(h-^ Mpc) 



ro'^ n dc 

(ft-' Mpc) (10-^ h^ Mpc-3) (/r' Mpc) 



0.25-0.75 0.53 136 -1.37±5]-^^ IS.TOzb^^' 1-00±046 l"'-40±3m 31.8±2.7 31.6±0.9 



0.75-1.50 0.97 141 -1.73±^ 



18.64±3 



1.21±" 



19.14±= 



!.96±0.76 48.1 ±1.4 



"Formal fits for rg and <5 were computed directly from the data. The error in A^ corresponds to the error in tq at the best-fit 5. 
''Best-fit parameters from frequentist fits. 
'^Best-fit ro from Bayesian marginalization. 



ULIRGs have space densities similar to the present cluster 
samples, a prerequisite for drawing evolutionary connections 
from these particular models. Thus the present work offers 
a measure of support for recent studies (e.g., Magliocchetti 
et al. 2007) indicating that ULIRGs may be associated with, 
or progenitors of, groups or low-mass clusters. 

5. CONCLUSIONS 

By deprojecting the angular correlation function measured 
in redshift bins spanning z = 0.25 to z = 1.5, we have deter- 
mined the real-space clustering amplitudes for ISCS clus- 
ters at (z) = 0.53 and (z) = 0.97 to be ro = 17.40!i];^^ and 
ro = 19.14;!;4 5^ /z~' Mpc, respectively. These measurements 
are consistent with the relation between correlation amplitude 
and mean intercluster distance predicted by LCDM. The ISCS 
clusters have total masses exceeding 10 '"^M© and will evolve 



into very massive clusters by the present day. 
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